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Abstract 

In this paper we prove the discrete compactness property for a wide 
class of p finite element approximations of non-elliptic variational eigen- 
value problems in two and three space dimensions. In a very general frame- 
work, we find sufficient conditions for the p- version of a generalized discrete 
compactness property, which is formulated in the setting of discrete differ- 
ential forms of order ^ on a polyhedral domain in M'^ (0 < £ < d). One of the 
main tools for the analysis is a recently introduced smoothed Poincare lift- 
ing operator [M. Costabel and A. Mcintosh, On Bogovskii and regularized 
Poincare integral operators for de Rham complexes on Lipschitz domains. 
Math. Z., (2010)]. In the case i = 1 our analysis shows that several widely 
used families of edge finite elements satisfy the discrete compactness prop- 
erty in p and hence provide convergent solutions to the Maxwell eigenvalue 
problem. In particular, Nedelec elements on triangles and tetrahedra (first 
and second kind) and on pai^allelograms and parallelepipeds (first kind) ai^e 
covered by our theory. 
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1 Introduction: Maxwell eigenvalue problem 



Maxwell's eigenvalue problem in a closed cavity G M'^ with perfectly con- 
ducting walls can be written as follows by means of the Maxwell-Ampere and 
Faraday laws: Find the resonance frequencies oj G M and the electromagnetic 
fields (E, H) ^ (0, 0) such that 

curl E = icj/iH 
curlH = —iueE 
E X n = 
H ■ n = 

where e and /i denote the dielectric permittivity and magnetic permeability, re- 
spectively. The fields E and H are sought in 

For simplicity, we consider now the case of homogeneous isotropic material 
with normalized material constants (e, /i = 1) — we will come back to the gen- 
eral setting in Remark 6.3. In a classical way, the elimination of the magnetic 
field from equations (1.1) yields the Maxwell eigenvalue problem with perfectly 
electrically conducting (PEC) walls in variational form: 

Seek u G H{cur\, fi) \ {0}, u; G Mg such that 

° (1 2) 

(curl u, curl v)^2(f^) = (u, v)^2(^) Vv G H{cur\,Q) . 

The elimination of the electric field would correspond to the same problem mod- 
elled through replacing H(cur\, Q) with i3"(curl, fi) ' . 

One aim of this paper is to prove the convergence of if (curl) -conforming 
Galerkin discretizations of Maxwell eigenvalue problem (1,2) in the framework 
of the p- version of the finite element method. The finite element approximation of 
Maxwell eigenvalues has been the object of intense investigations for more than 
20 years. It was soon recognized that the if(curl)-conforming Galerkin finite el- 
ement discretizations need special finite element spaces that are generally termed 
edge finite elements (see [44, 45, 14]). 

The first attempts to analyze the discretized eigenvalue problem have been 
made for the /i-version of edge finite elements. We mention [39] as a pioneer- 
ing work on lowest order edge finite elements, where the discrete compactness 
property (see [2]) has been indicated as a key ingredient for the analysis. Other 
relevant works on the subject are [13, 8, 19, 43, 40, 24, 9], and we refer the inter- 
ested reader to [37, 42] and to the references therein for a review on this topic. 

In these references, the Maxwell eigenvalue problem is often studied using 
variational formulations different from (1.2), for example mixed formulations [9], 

'By and large, we adopt the standard notations for Sobolev spaces, see [34, Ch. 2]. 



in il 
in VL 
on dVL 
on dVt, 



(1.1) 
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regularized formulations [23, 25] or mixed regularized formulations [5, 17]. With 
the exception of the method of weighted regularization [23, 25, 17], where H^- 
conforming elements can be used, these formulations use the i7(curl)-conforming 
edge elements. In their analysis, special conditions implying convergence of the 
discrete eigenvalue problems are presented, for example the so-called Fortid prop- 
erty [8], or the GAP property [16]. As explained there, these conditions are related 
to the discrete compactness property. Here we choose to work with the simple 
variational formulation (1.2) and its generalization to differential forms. The role 
of the discrete compactness property in this context has been discussed in detail 
in [19]. 

The analysis presented in the references above covers the /i-version for ba- 
sically all known families of edge finite elements. It soon turned out, however, 
that the analysis of the p- and /ip-versions of edge finite elements needed tools 
different from those developed for the /i-version. In [12] the two-dimensional tri- 
angular case has been studied for the hp-version, but the analysis depends on a 
conjectured estimate which has only been demonstrated numerically. In [II] a 
rigorous proof for the /ip- version of 2D rectangular edge elements has been pro- 
posed (allowing for one-irregular hanging nodes) which, in particular, contains the 
first proof of eigenvalue/eigenfunction convergence for the pure spectral method 
(p- version with one element) on a rectangle. 

What paved the way for a successful attack on a general p-version analysis 
was the regularized Poincare lifting recently introduced in [26]: it enjoys excellent 
continuity properties and at the same time respects discrete differential forms. In 
this paper we are going to show how the regularized Poincare lifting can be com- 
bined with another recent invention, the projection based interpolation operators, 
see [27, 29], to clinch the analysis of the p-version of edge elements. This al- 
lows to prove the discrete compactness (and hence the convergence of the discrete 
eigensolutions) for a wide class of finite elements related to discrete differential 
forms: for (1.2) this includes, in particular, Nedelec elements on triangles and 
tetrahedra (first and second kind) and on parallelograms and parallelepipeds (first 
kind). 

As already mentioned, one of the key ingredients for the convergence analysis 
is the discrete compactness property. Much insight can be gained from investi- 
gating it in the more general framework of discrete differential forms (see [4] for 
a lucid introduction to this subject). In this setting, the proofs are more natu- 
ral and simultaneously cover, in particular, two- and three-dimensional Maxwell 
eigenvalue problems. 

Plan of the paper. The structure of the paper is as follows. We start in Sec- 
tion 2 with a generalization of (1.2) to eigenvalue problems associated with the 
de Rham complex on differential forms. Then we define the discrete compactness 
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property and discuss its significance in the context of Galerkin discretization: in 
association with two standard completeness properties, it gives a crucial sufficient 
condition for the convergence of eigenvalues and eigenvectors. Section 3 is the 
core of our paper and contains the description of our abstract assumptions. Hav- 
ing in mind the p- version of finite elements, we consider a fixed mesh of a 
bounded Lipschitz polyhedron C M'' and a sequence of spaces of discrete dif- 
ferential forms of order i (with < i < d) together with projection operators 
onto discrete spaces; we prove that our assumptions imply the validity of the dis- 
crete compactness property for such a sequence of spaces (Theorem 3.2). The 
abstract theory relies on the existence of suitable Poincare lifting operators which 
are presented in Section 4. The mapping properties of these lifting operators allow 
to specify some of the function spaces appearing in our abstract assumptions. In 
Section 5 we recall the classical families of discrete differential forms with high 
degree polynomial coefficients on simplicial or tensor product elements. 

Our abstract theory applies to any dimension d, but for want of suitable regu- 
larity results, embeddings, and projection operators, we can give examples satis- 
fying all of its assumptions only in dimensions d = 2 and d = 3. This is done in 
Section 6, where we concretize the function spaces and recall embedding results 
and properties of projection based interpolation operators related to these spaces. 
All abstract assumptions are then satisfied, leading to the main convergence result 
stated in Theorem 6. 1. The analysis of a p- version edge element discretization of 
the Maxwell eigenvalue problem (1.2) is covered as case d = 3 and i = 1, see 
Corollary 6.2. 

2 Differential forms and generalized Maxwell eigen- 
value problem 

The variational eigenvalue problem (1.2) turns out to be a member of a larger 
family of eigenvalue problems, when viewed from the perspective of differen- 
tial forms. This more general perspective offers the benefit of a unified theoretical 
treatment of different kinds of eigenvalue problems, e.g., the scalar Laplace eigen- 
problem. Maxwell cavity eigenproblems in dimensions 2 and 3, the eigenproblem 
for the grad div-operator in dimension 3. This policy has had remarkable success 
in numerical analysis recently, cf. [3]. Thus, in this section we first recall some 
basic notions related to differential forms. We refer the interested reader to [4, 
Sect. 2] for an introduction to this subject. 
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2.1 Function spaces of differential forms 

Given a bounded Lipschitz domain Q C M°', we denote by C°°{^1, A^), < i < d, 
the space of smooth differential forms on and by A(f2) the corresponding anti- 
commutative graded algebra 

e 

The exterior derivative is a graded linear operator of degree one 

d : A(fi) ^ A(f]), 

that is, for any £ it is represented by : A^ ^ A^+i). 

We rely on the Hilbert spaces 

H{di,n) ■= {v G A^ : d^v G /.^(l], A^+i)} , (2.1) 

where L^(n, A^) is the space of differential ^-forms on with square integrable 
coefficients in their canonical basis representation, see [26, Sect. 2]. Its inner 
product can be expressed as 

(u,v)o^:= / uA^v, u,YeL\Q,A'), (2.2) 
Jn 

with ^ the Hodge star operator induced by the Euclidean metric on R'^, which 
maps ^- forms to (d — £) -forms. As above, a o tags the subspaces of forms with 
vanishing trace tiQQ on d^l, which can also be obtained by the completion of 
compactly supported smooth f-forms with respect to the H^d^, r2)-norm: 

if (d,, n) := {v G Hid,, n) : tian v = 0}. (2.3) 

o 

The subspace of closed forms is the kernel of de and is denoted by H{de 0, Q): 

H{de 0, n) := {v G H{de, (]) : d^ v = 0}. (2.4) 

Finally, we adopt the standard notation for the exterior co-derivative operator 
6e :=^dd-£^; in particular, we have 6e : C°°(fi,A^) ^ C°°(fi,A^-i). 

2.2 Variational eigenvalue problems 

After choosing bases for the spaces of alternating multilinear forms on MJ^, vector 

fields ("vector proxies") Q i— )■ M^*^) provide an isomorphic model for differential i- 
forms on il. Choosing the standard "Euclidean basis", the operators -k, 5, tTgn are 
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Table 1: Identification between (operators on) differential forms and (operators 
on) Euclidean vector proxies in and M'^ 



Differential form 


Proxy representation 




d = 2 


d = 3 


£ = 


do 


grad 

<P\dn 


grad 


1=1 


di 

Si 


curl 

(u X n)|an 
H{cur\,n) 
div 


curl 

(u X n)|an 

H{cur\,n) 

div 


i = 2 


d2 

52 





curl 


div 

(q- n)|oc 

H{6\v,n) 

curl 



incarnated by familiar operators of classical vector analysis, different for different 
dimension d and degree i, see Table 1 and [4, Table 2.1]. 

Hence, the eigenvalue problem (1.2) with e, /i = 1 is the special case d = 3, 
i = 1, of the following variational eigenvalue problem for differential £-forms, 
0<i<d: 

Seek u e H{de, fi) \ {0}, lo e M^, such that 

° (2 5) 

(df u, d^ v)o = (u, v)g W e H{de,n) . 

A key observation is that the bilinear form (u, v) i— (d^ u, d^ v)^ ^ has an infi- 

o 

nite dimensional kernel H{diO,Q) comprising all closed ^-forms. It provides the 
invariant subspace associated with the essential spectrum {0} of (2.5). This es- 
sential spectrum can be identified as the main source of difficulties confronted in 
the Galerkin discretization of (2.5). 

On the other hand, any solution u of (2.5) for a; 7^ satisfies (u, d^^i i/>)q ^ = 

o 

for all '0 G if (d£_i 0, Q). Thus the eigenfunctions corresponding to non-zero 
eigenvalues belong to the subspace 

Y{de, n) := {v e H{de, Q) : (v, d,_i V)o,n = V^^ G ^(d,„i, Q)}, (2.6) 
which means they belong to the kernel of 5i. This is the generalization of the 
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divergence free constraint found for electric fields in the Maxwell case. From 
[46] we learn the following theorem. 

Theorem 2.1 For any d E N, < I < d, the embedding 
is compact. 

o 

By restricting the eigenvalue problem to F(d^, we can therefore use Riesz- 
Schauder theory. This implies that (2.5) gives rise to an unbounded sequence of 
positive eigenvalues = {u'^Y 

AO = < < < . . . , A*^ ^ oo (A; ^ oo) , (2.7) 

with associated finite dimensional mutually L^(fi) -orthogonal eigenspaces. 

Remark 2.2 Owing to the zero trace boundary conditions imposed on the func- 
tions in (2.5), it may be called a Dirichlet eigenvalue problem. Using H{di, f2) 
as variational space would result in the corresponding Neumann eigenvalue prob- 
lem. Its analysis runs utterly parallel to the Dirichlet case using the techniques 
presented below. □ 



2.3 Approximation of the eigenvalue problem and the role of 
discrete compactness 

In the sequel we fix the degree i, < i < d, of the differential forms. Spaces of 
discrete differential forms 

Vl C H{di, Q) , dim < oo , 

lend themselves to a straightforward discretization of (2.5). In this section, p eN 
stands for an abstract discretization parameter, and, sloppily speaking, large values 
of p hint at trial/test spaces of high resolution. 

We consider the approximation of the eigenvalue problem (2.5) by the Galerkin 
method: 

Find UpEVi\ {0}, uj E M([, such that 

o (2 8) 

(d^ Up, d£ Vp)g_^ = (up, Vp)g_^ Vv G Vj . 

Now, the key issue is convergence of eigenvalues and eigenvectors as p — )• oo, 
rigorously cast into the concept of spectrally correct, spurious-free approximation 
[19, Sect. 4]. Let us recall these notions in a few words for the case of self-adjoint 
nonnegative operators without continuous spectrum (which is the case here). 
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The spectral correctness of the approximation of an eigenvalue problem such 
as (2.5) by a sequence of finite rank eigenvalue problems (2.8) means that all 
eigenvalues and all eigenvectors of (2.5) are approached by the eigenvalues and 
eigenvectors of (2.8) as p — oo. If (2.5) has a compact resolvent (which is the 
case only when £ = 0), the spectral correctness is an optimal notion: It implies 
that if {A'^}fc>i and {\p}k>i are the increasing eigenvalue sequences of (2.5) and 
(2.8) (with eigenvalues repeated according to their multiplicities), then 

Aj ^ A'^ as p-^oo \fk>l, (2.9) 

and the gaps between eigenspaces (correctly assembled according to multiplicities 
of the eigenvalues of (2.5)) tend to as p — j- oo. 

If we face an eigenvalue problem for a self-adjoint non-negative operator with 
an infinite dimensional kernel, and otherwise discrete positive spectrum (which is 
the case for (2.5) for all £ > 1), the spectral correctness implies the same prop- 
erties as above with the following modifications of the definitions: Now {X''}k>i 
is the increasing sequence of positive eigenvalues of (2.5) (as specified in (2.7)) 
and, given a positive number e < A\ {Ap}fc>i is the increasing sequence of the 
eigenvalues of (2.8) larger than e (still with repetitions according to multiplici- 
ties). With such conventions, spectral correctness still implies convergence of 
eigenvalues (2.9) and eigenspaces as above. In this context, spurious-free approx- 
imation means that there exists 5o > such that all eigenvalues of (2.8) less than 
Eq are zero. Therefore, spectrally correct, spurious-free approximation implies the 
convergence property (2.9) and the corresponding convergence of eigenspaces, if 
we define {\p}k>i as the increasing sequence of the positive eigenvalues of (2.8). 

There exist several different ways, all well studied and summarized in the lit- 
erature of the last decade, for proving the convergence of the discrete eigenvalue 
problem (2.8) to the continuous eigenvalue problem (2.5): One can use a reformu- 
lation as an eigenvalue problem in mixed form as analyzed in [9], or one can use a 
regularization which gives an elliptic eigenvalue problem for the Hodge-Laplace 
operator as analyzed in [4], or one can follow the arguments of [ 1 9] and study the 
non-elliptic problem (2.5) directly. 

Here we outline the latter approach, which employs the analysis of [32] of 
the approximation of eigenvalue problems of non-compact selfadjoint operators. 
Since [19] deals only with the Maxwell case, i. e. d = 3, £ = 1, we examine 
the main arguments, in order to verify that they are also valid for the general case. 
The proofs we give are adaptations of those of [19] to our more general situation. 

Let us define the solution operator A : L^(f2, A^) — )■ i^(d^, Vt) of the source 
problem corresponding to the eigenvalue problem (2.5) and its discrete counter- 
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part Ap-.V^-^ by 

(d, Ai, d, v)o + (Af , v)o = (f , v)o,^ Vv G #(d,, fi) 
(dMpf,d,v)o^^ + (A,f,v)o^^ = (f,v)o^^ VvGVj. 

Note that the operators A and have the same eigenfunctions and the same 
eigenvalues (after a transformation) as the eigenvalue problems (2.5) and (2.8). 
Namely, (2.5) and (2.8) are equivalent to the relations 

u = (w^ + l)ylu ; Up = {uj'^ + l)ApUp. (2.11) 

The infinite-dimensional eigenspace at u; = shows that A is not a compact 
operator. 

Following [19], three conditions are identified that together are necessary and 
sufficient for a spectrally correct, spurious-free approximation of A by Ap or, 
equivalently, of the eigenvalue problem (2.5) by the discrete eigenvalue problem 
(2.8). 

The first condition is rather natural. It states that the sequence of discrete 
spaces (Vp)pgj^ is asymptotically dense in H{de,il) (compare [19, Condition 
(CAS) - completeness of approximating subspaces]) 

(CAS) lim inf ||v- Vp||j:^(j^_^) = W E H{di,n) . (2.12) 

The second condition, only relevant for £ > 0, states that closed forms can be 
well approximated by discrete closed forms (compare [19, Condition (CDK) - 
completeness of discrete kernels]) 

(CDK) lim inf ||z - Zp||^2(m = Vz G iJ(d^ 0, 1]) . (2.13) 

The third condition is the most intricate one and has been dubbed discrete com- 
pactness. For its formulation, we introduce the orthogonal complement space of 
the discrete closed forms: 

ij:={upGVj: (up,Zp)o^^ = Vzp G n if (d, 0, fi)}. (2.14) 

Definition 2.3 Let us choose £ G {1, . . . , c? — 1}. The discrete compactness prop- 
erty holds for a family (Vp )^^^ of finite dimensional subspaces of H{di, f2), if for 
any subsequence N' o/N, any bounded sequence 

(up)pgj^, C H{di, n) with Up G Zl 
contains a subsequence that converges in L'^{Q, A^). 
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The convergence proof is based on two lemmas, the first of which corresponds 
to [19, Theorem 4.12]. It implies, according to [32, Condition PI) and Theorems 
2,4,5,6], the spectral correctness of the approximation. 

Lemma 2.4 If (2.12) and the discrete compactness property hold, then 

lim sup \\A\p - ^pVpll^. , m = • (2.15) 

VpeV^;||vp||j:^(^^_s^)=l 

Proof. Note first that for Vp G fl H{6i 0, Vt) there holds Avp = Vp = ApWp, 

o . 

so that by orthogonal decomposition of one gets 

sup W^^p- ^p^p\\H(d,,n) = ^ sup pVp-ApVp||^(j^^^) . 

Furthermore, one has by definition of A and Ap 

\\Awp - ApVp||^(,^^^) = inf \\Awp - Wp||^(,^^^) . 

Assume now that (2.15) does not hold. Then there exists e > 0, a subsequence M' 
of N and a sequence (vp)pgN' with Vp G Zp satisfying ||vp||^^j^ n) ^ -'- ^'^^ 

Pvp - wp||^(,^^^^) >e Vp G N', wp G . (2.16) 

We can apply the discrete compactness property to the sequence (vp) and ob- 
tain a subsequence converging in A^') to some v G A^). Since A : 
L^(r2,A^) — )■ H(d£,il) is continuous, we find G H{de,^l), and the approx- 

o ^ 

imation property (2.12) provides us with a sequence (wp) with Wp G Vp that 

O 

converges in H{di, Q) to Ay. Hence for the subsequence we obtain 

PVp - Wp||^(,^^^) < PVp - Av||^(,^ + ||AV - Wp||^(,^^^) ^ , 

in contradiction with (2. 1 6). □ 
The second lemma corresponds to [19, Corollary 2.20]. It gives the discrete 
Friedrichs inequality (in [9] also called "ellipticity in the discrete kernel"), and it 
is easy to see that this implies that w = is not a limit point of positive discrete 
eigenvalues, so that the spurious-free property of the approximation follows. 

Lemma 2.5 If (2.13) and the discrete compactness property hold, then there ex- 
ists a > such that for all p eN 

Vv G 4 (2.17) 



10 



Proof. Assume that (2.17) does not hold. Then there exists a subsequence N' of 
N and a sequence (vp)pgN/ with Vp G Z^^ satisfying 

l|vpllL2(f^) = 1 and ^lim ||d^Vp||^2(^) = 0. (2.18) 

The discrete compactness property can be applied to this sequence and gives a sub- 
sequence converging in L^(f2, A^) to some z G A^). From (2.18) follows 

o o 

that the convergence actually takes place in H{6i,Vt) and that z G iJ(d£0,fi). 
Therefore the approximation property (2.13) provides us with a sequence (zp) 
with Zp G Vp n H{6i 0, Vt) that converges in LF'{Vt, A^) to z. Hence for the subse- 
quence we find 

o „ o . o 

This leads to a contradiction, because Vp G Zp and Zp G fl if(d^0,r2) are 
L^(r2) -orthogonal, hence for all p 

W'^P ~ ^pllL2(n) = ll'^pllL2(n) + NpllL2(n) > 1 • 

□ 

To summarize. Lemmas 2.4 and 2.5 together prove the following result. 

Theorem 2.6 If the completeness of approximating subspaces (2.12), the com- 
pleteness of discrete kernels (2.13) and the discrete compactness property hold, 
then (2.8) provides a spectrally correct, spurious-free approximation of the eigen- 
value problem (2.5). 

Remark 2.7 The main focus of this section is on the convergence of the eigen- 
values and the eigenfunctions of problem (2.8) to those of (2.5). On the other 
hand, when considering concrete applications it is crucial to investigate the order 
of convergence. In order to do so, several strategies are available. A straightfor- 
ward approach which well fits the theory summarized in this section makes use 
of the results from [33]. Theorem 1 of [33] states in this particular situation that 
the error in the eigenfunctions (measured as usual by the gap of Hilbert spaces) 
is bounded by the best approximation, and Theorem 3(c) of [33] states that the 
eigenvalues achieve double order of convergence since our problem is symmetric. 
An alternative approach makes use of the equivalence of problems (2.5) and (2.8) 
with suitable mixed formulations [10, Part 4]; in this case an estimate of the or- 
der of convergence can be achieved by the standard Babuska-0 shorn theory for 
the spectral approximation of compact operators applied to the mixed formula- 
tions [10, Theorems 13.8, 13.10, 14.9, 14.11]. □ 



11 



3 An abstract framework implying discrete compact- 
ness 



In this section we fix a degree of differential forms 

£e{l,...,d-l}, 

and we formulate a set of hypotheses which allow us to prove the discrete com- 
pactness property. These hypotheses are organized in three groups: 

o . 

1. standard assumptions related to the finite element spaces (Sect. 3.1), 

2. assumptions on the existence and key properties of "lifting operators" (Sect. 3.3), 

o , 

3. hypotheses on projections onto complying with the commuting diagram 
property and satisfying an approximation property (Sect. 3.4). 

To state these assumptions we have to introduce intermediate spaces X and S of 
more regular forms 

Vl C X{m, AO C H{de, n) and Vj-^ c S{m, A^^^ C i^(d^_i, Q) , 

allowing compact embedding arguments and precise notions of continuity of lift- 
ing and projection operators. 

3.1 Discrete spaces 

Our focus is on finite element spaces. For the sake of simplicity, we restrict our- 
selves to polyhedral Lipschitz domains fi. We assume that the finite dimensional 
trial and test spaces V^, p e N, are based on a. fixed finite partition DJl of fi, 
composed of elements (cells) K: 

n=[jK , K nK' = (ll , if K ^ K', K,K' em . 
Kern 

For a cell K E VJl, let ^rn.(,K) designate the set of m-dimensional facets of K: 
for m = these are the vertices, for m = 1 the edges, for m = d — 1 the faces, 
mdUK) = {K}. 

We take for granted that the discrete spaces Vp can be assembled from local 
contributions in the sense that for each mesh cell K em there is a space Vp{K) C 
C^{K, A^) of smooth forms on K, such that 

Vj = Vim := {v G if(d,,fi) : v|^ G V,^(/i) 'iK em] . (3.1) 
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In other words, can be defined by specifying the local spaces Vp{K) and re- 
quiring the continuity of traces across inter-element boundaries as well as the 
boundary conditions on d^l. 

In the same fashion, we introduce a corresponding family Vp~^ C iJ(d^_i, f2) 
of spaces of discrete (£ — 1) -forms. We will see later on that as a consequence of 
further hypotheses, the local spaces Vp~^ (K) and Vp{K) satisfy an exact sequence 
property. 

3.2 Spaces of more regular forms 

We introduce a Hilbert space X{Vyt, A^) C H{di,Q) that captures the extra regu- 

O 

larity that distinguishes f-forms in the space Y{di, Q). We can think of this space 
as a space of "more regular" £-forms on 



Assumption 1 The space F(d£, fl) defined in (2.6) is continuously embedded in 

x{m, My. 



This means that with C > depending only on Q 



X(OT,Af) 



<C||u||^(,^^^) VuGr(d,,fi). (3.2) 



On the other hand, X{Dyt, A^) has to be small enough to maintain the compact 

o 

embedding satisfied by F(d^, Q), cf. Thm. 2.1. 



Assumption 2 The space X{Tl, M) is compactly embedded in A^) 

X(9Jt,A0 'TL2(fi,A0. 



As with the discrete spaces, the spaces X(9Jt, A^) are built from local contri- 
butions and will therefore depend on the mesh 97t. We assume that for each mesh 
cell K E'^XR there are Hilbert spaces X(K, A^) so that: 

X{m, AO = {v G H{de, n): v |^ G X{K, A^ VA' G Tl} , (3.3) 
and, in addition, the norm of X{Tl, M) is defined through local contributions: 

1 1 2 1 1 2 ^ — ^ II I 1 1 2 



Finally, the local spaces have to be large enough to contain the discrete forms for 
any value of p: 

V'p{K)cX{K,M) . (3.5) 
13 



Correspondingly, we introduce a space 5(971, A^~^) C i7(df„i,r2) of "more 
regular potentials". Similar to X{Tl, A^), the spaces S{9Jl, A^~^) are mesh-depen- 
dent and allow for a characterization through local Hilbert spaces ^(A', A^~^), 

K em, 

S{m, A^-i) = {-0 G H{de^i, Q): tl^lj^e S{K, A^-i) eM] . (3.6) 
They are endowed with the norm 

ll^llW^-) = ll^llWi.f^) + Yl ll'^lxll5(x,A^-i) • (3.7) 

The local spaces are large enough to contain the local discrete potential spaces: 

V'p-\K) c S{K,A'-^) . (3.8) 
The following assumption establishes the connection between X{'[Dl, A^) and 

s{m, A^-i). 

Assumptions The exterior derivative maps 5(971, A^"^) continuously into 

x{m, A^y. 

and the image is maximal: 

d^_i S{m, A^"^) = d^„i H{de^i, Q) n x(mt, A^). 



To conclude this subsection, note that in the case of an element K touching the 
boundary dil, like for the discrete spaces Vp{K) and Vp~^{K), the local spaces 
X{K,A^) and S{K,A^~^) are not obliged to comply with any boundary condi- 
tions. 

3.3 Local liftings 

A pair of linear mappings Rk,K ■ C°°{K,A'') ^ C~(A', A^-^), k = i,£ + l,is 
called a lifting operator of degree I if it fulfills 

d^_i o R^^^ + R£+i,i^ o d^ = ld£ . (3.9) 

This relation characterizes a "contracting homotopy" of the de Rham complex [5, 
Section 5.1.2]. 

Besides this algebraic relationship, our approach hinges on smoothing proper- 
ties of the lifting operators, expressed by means of the local spaces S{K, A^~^) of 
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more regular potentials and X{K^ A^) of more regular forms. The next assump- 
tion summarizes the continuity expected from the lifting operator. 

Assumption 4 For every A' G SOT there is a lifting operator (R^^x, ^(.+i,k) whose 
components can be extended to continuous mappings 

R,+i,j^ : L\K, A^+i) ^ X{K, and R^,k : X{K, K') ^ S{K, A'-') , 

and thus identity (3.9) holds on X{K, A^). 

As a consequence, for each cell K G 971, we have the exact sequence 

S{K,A'~^) X{K,A') L\K,A'+^). (3-10) 

Finally, the local liftings have to be compatible with the local spaces of discrete 
differential forms: 

Assumption 5 The local operators R^+i z^, when applied to exact local discrete 
{£ + l)-forms, yield local discrete ^-forms, i.e.. 



3.4 Local projectors 

As usual in methods based on discrete commuting diagrams we need projection 
operators tt^ j^- onto discrete spaces for (i — 1) -forms and forms. For degree 
i — 1, our local spaces S{K, A^~^) of more regular potentials can play the role of 
domains for the projectors vr^"^. For the degree i, by generalization of what we 
actually need in the case of dimension d = 2 and d = 3 for Maxwell, we define 
our projectors vTp on smaller spaces than X(A', A^). We denote these new spaces 
by S{K, A^) and require that they contain for all p the p-dependent subspaces 

Xp{K,A') = {ueX{K,A'):deuedeV'piK)}. (3.11) 

On the same model as (3.6)-(3.7), we define the corresponding global spaces 
S{m, M) and 

Xp{m,A^) = {u G X{m,A^) : d^u G deV^p} (3.12) 
and we have the continuous embeddings 

Xp{Dn, A^) ^ 5(531, A^) ^ X{m, A^) . (3.13) 
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Assumption 6 There are local continuous linear projections 

vrj;^ : S{K, A'-') ^ Vf^K) and <^ : S{K, A') ^ V^K) 
for all mesh cells K E VJt. 

The standard commuting diagram property is as follows. 

Assumption 7 The projectors tt^^ and n^ j^ are compatible with the exterior 
derivative in the sense that the diagram 

S{K,A'-') S{K,A') 

commutes for every K E OJt. 

Let us note that, as a consequence of Assumptions 4 and 7, we find that the 
sequence 

VtHK) ^ V^(ir) d, (V^(JO) 

is exact. 

Besides, the local projections acting on (i — l)-forms are supposed to enjoy a 
crucial approximation property using the Hilbert space norms || ■ H^^^. a^-i)- 

Assumption 8 There is a function ei^i : N (-t- R'^ with lim ei^i{p) = so that 

p— >oo 

Finally we assume for the projections tt^ ^ a natural condition of conformity: 
For all u E Xp{K, 

trj.u = => tTF7rlKU = E dm{K) , i < m < d , (3.14) 

and the corresponding condition for the projections tt^"^. This makes it possible 
to define global linear projections 

Til: S{m,A')^Vl and vrj"^ : 5(971, A^ ^ Vj"' 

by patching together the local operators 

(7r^u)|^:=<^(u|j and « V) |^ := vr^7;,i(0 |j V/TgOH. (3.15) 
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As a consequence of Assumption 7 and (3.15), the global projectors tt^ ^ and vr^ 
inherit the global commuting diagram property 



s{m,A'-') Sim, A') 

i (3.16) 



3.5 Proof of the discrete compactness property 

The estimate of Assumption 8 on "potentials" carries over to ^-forms with a dis- 
crete exterior derivative, that is, the elements of the space Xp{dR, A^), see (3.12). 

Lemma 3.1 (Global projection error estimate) Making Assumptions 4 through 8, 
the estimate 

\\^~^t^\\LHn,M) - C^i-i(p) l|u|lx(OT,AO Vu e Xp(9Jl,A^) 

holds true, with a constant C > independent of p. 

Proof. Pick any u e Xp(97l, A^). The locality of the projector vr^, cf. (3.15), 
and (3.4) allow purely local considerations. Single out one cell K E Tl, still write 
u = u 1^ G Xp(K, A^), and split u on K using (3.9) from Assumption 4: 

u = d^_i Rej^u + R£+i,/< d£ u = d^_i + Re+i,K de u . (3.17) 

with <p := i^u. The continuity of Ri^k from Assumption 4 reveals that 

U\\s{K,A'-') ^'^Mx{KA') ' (3-18) 

where here and below C will denote constants (possibly different at different oc- 
currences) which depend neither on u nor on p. 

Thanks to identity (3.17) and the commuting diagram property from Assump- 
tion 7, we have 

■Kp j^u = d£_i 7ip-jl^(f) + 7ip j^Ri+i,K di u . (3.19) 

Recall that u G Xp(K, A^) belongs to the domain of vr^^ by Assumption 6. Fur- 
ther, as u G Xp{K, A^), from Assumption 5 we infer that 

R,+^,KdiueV'p{K) . (3.20) 

Thus, owing to the identities (3.17), (3.19) and the projector property of tt^^, the 
task is reduced to an interpolation estimate for iTp'j^: 

(Id - n'p j,)u = d,_i(ld - 7rJ-^)0 + (Id - irl ;,)R^+i,x d, u . (3.21) 

^ ' 

=0 by (3.20) 

17 



As a consequence, invoking Assumption 8, 



(3.18) 

<ee-i{p)\\4)\\s(K,Ai-^) < C'^^-ib) l|u|lx(i^,A*) ' (3-22) 

which furnishes a local version of the estimate. This estimate is uniform in K E 
Tl because 9Jt is finite. Due to (3.4), squaring (3.22) and summing over all cells 
finishes the proof. □ 
We are now in the position to prove the main result of this section. 

Theorem 3.2 (Discrete compactness) Under Assumptions 1 through 8, the dis- 
crete compactness property of Definition 2.3 holds for the family ( of sub- 

O 

spaces of H{di, Q). 

Proof. The proof resorts to the "standard policy" for tackling the problem of 
discrete compactness, introduced by Kikuchi [39, 40] for analyzing the /i-version 
of Whitney- 1 -forms. It forms the core of most papers considering the issue of dis- 
crete compactness, see [12, Thm. 2], [11, Thm. 11], [37, Thm. 4.9], [31, Thm. 2], 
etc. 

O 

Let us introduce the discrete analogue of the space ^(d^, 

Z ■■= {vp e Vl : (v„ d,_i Vp)o,f, = VV^, G V'-'}. (3.23) 

The space contains as a subspace. 

We consider a subsequence N' of N and a H{di, f2)-bounded sequence (up)^^^, 

with members in Z^. Thus Up belongs in particular to 3^* and the sequence 
("p)peN' satisfies 

(i) UpeV'p, (3.24) 

(ii) (up,d,_iVp)o,f, = VV'pGVJ"^ (3.25) 

(iii) l|upL(d„n)<l VpeN'. (3.26) 

We have to confirm that it possesses a subsequence that converges in A^). 
We start with the L^(f2, A^)-orthogonal projection of Up into Y{di,^l) and 

o o 

parallel to d^-i H{di^i, Q): let Up be the unique vector field in H{de, ft) with 

Up = Up + d^„i 0p, E H{d(,^i, n) , (3.27) 
(Gp, d,_i V')o n = G i7(d,_i, n) . (3.28) 
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Obviously, the latter condition implies 

Up e Y{6e, n) . (3.29) 

Hence, by virtue of Assumption 1, the fact that d^Up = d^Up, and (3.12), Up 
satisfies 

Up e Xp{n,A ), ||up||^(2jj^^£-) < ||up||^(j^ , (3.30) 

where C > does not depend on p. 

Since d^_i (/)p = Up — Up E X{Tl, A^), Assumption 3 implies that we may 

assume that 4>p G S{m, A^^^). 

Thus we can use Nedelec's trick [44] to obtain 

11^ Il2 / — I — I'-' \ 

W^P ~ %IIl2(q^a«) ^ l^P ~ Up7Up — TTpUp + VTpUp — Up) 

= [Up-Up,Up-npUp)^^^ . 
This holds because from (3.27) and the projector property of tt^ we know 

TTpUp - Up = VTpUp + VTp d^_i 4>p- Up = VT^ di-i (f)p , 

and thanks to the commuting diagram property (3.16) (deduced from Assump- 
tion 7) combined with the orthogonality conditions (3.25) and (3.28), we find 

(up - Up, TTpUp - Up)^^^ = (up - Up, de-i vTp-^^pj^ ^ = . (3.32) 
Hence, appealing to Lemma 3.1, with C > independent of p, we get 

II^P ~ ^pllL2(n,A«) — II^P ~ ^P^plLsfQ A*^) — ll"pllx(OT,A«) 

(3.3„) ' (3.33) 

From (3.30) we conclude that the sequence (up)^^^^, is uniformly bounded in 
X(9Jl, A^). By Assumption 2 it has a convergent subsequence in L^(r2, A^). Ow- 
ing to (3.33), the same subsequence of (up)pgj^, will converge in L^(f2, A^). □ 

3.6 Approximation of the eigenvalue problem 

As discussed in Section 2.3, the discrete compactness property is the cornerstone 
of the proof of the convergence of the discrete generalized Maxwell eigenvalue 
problem (2.8). 
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Corollary 3.3 In addition to the hypotheses of Theorem 3.2, namely Assump- 
tions 1 through 8, assume that property (CAS) (2.12) holds and that the space 
X(Wl, A ) n H{di 0, Q) is dense in H{di 0, Q). Then (2.8) provides a spectrally 
correct, spurious-free approximation of the eigenvalue problem (2.5). 

Proof. We use Theorem 2.6 from Section 2.3. Considering that the dis- 
crete compactness property is provided by Theorem 3.2, and that we assume the 
approximation property (CAS) (2.12), we only need to show the approximation 
property (CDK) (2.13), which concerns the approximation of closed forms by 
closed discrete forms. 

Since we assumed the density of X{Wl, A^) fl H{de 0, Q) in H{dt 0, Q), it is 
sufficient to prove (CDK) for z e X(9Jt,A^) fl H{diO,n). Such z belongs to 
XpidJl, A^), and we can therefore apply Lemma 3.1, which shows that vr^z — )■ z 
in A^). We will have accomplished to show (CDK) with Zp = vr^z, as soon 

as we show that Zp = 0. Keeping in mind that Zp e Vp C H{di, Q), we see that 
it is sufficient to show the local relation d^ Zp = in K for every K G Tl. This 
follows finally as in (3 . 1 9) in the proof of Lemma 3.1, because d^ z = implies 

Hence Zp = de tt^ ^^z = d^ d£_i TTp~^Ri^K'^ = 0, which ends the proof. □ 

Remark 3.4 The abstract theory developed in this section can be applied to the 
/i- version of discrete differential forms, if the dependence of the constants on the 
size of the cell K is made explicit by means of scaling arguments. Here, we forgo 
this extra technicality and refer the reader to [37, Sect. 4.4]. □ 

4 Regularized Poincare lifting 

In this section we describe the construction of a local lifting operator that will 
satisfy Assumptions 4 and 5 in Section 3.3 for suitable spaces X{K, A^), S{K, A^) 
and Vp{K). We follow the presentation in [26], where these operators are analyzed 
and where it is shown in particular that they are pseudodifferential operators of 
order —1. 

4.1 Definition 

We consider a bounded domain D C that is star-shaped with respect to some 
subdomain B C D, that is, 

MaeB.xeD: {(1 - t)a + tx, < t < 1} C D . (4.1) 
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For a G B and 1 < ^ < rf, we define the Poincare operator acting on a 
differential form u G C°°{D, A^), by the path integral 

R^,au(x) = (x - a) J /" t*^^^ u{a + t{x - a)) dt , xeD. (4.2) 
Jo 

Here the symbol j denotes the contraction (also called "interior product") of the 
vector field x ^-^ (x — a) with the £-form u. It is clear that R^ « maps C°°{D, A^) 
to C'^{D, A^~^) and it has been shown (see [35] for proofs in the case d = 3) that 
it can be extended to a bounded operator from L'^{D, A^) to L'^{D, A^"^). In order 
to define the regularized Poincare operator R^, we choose a function 

e E Co°^(M"') , supp 9c B, [ e{a) da = I , 

Jb 

and set 

Reu(x) = / 9{a)Re^au{x) da . (4.3) 
Jb 



4.2 Regularity 

The substitution y = a + t(x — a), t = — t) transforms the double integral 
in (4.2), (4.3) into 



oo 



R^u(x) = j j {t- ly-^T'^-^e {x + T{y - x)) (x - y) J u{y) dr dy 

1 (4.4) 

where the kernel k{y, z) has an expansion into quasi-homogeneous terms: 



/ s'-\s + lY-'e{y + sz)ds 
Jo 



k{y^z) = -z 

j=o \ \ Jo II 

The operator R^ is therefore a weakly singular integral operator. In [26, Section 
3.3], the following result is shown. 

Proposition 4.1 For 1 < i < d, the operator Ri is a pseudodijferential oper- 
ator of order -1 on W^. It is well defined on C°°{D,A^), it maps C°°{D,A^) 
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to C^iD, A^-i) and C°^{D, to C°°{D, A^-^), and for any s e ^ it has an 
extension as a bounded operator 

Here, H^{D, A^) is the Sobolev space ofi-forms on D of order s. 
4.3 Lifting property 

The lifting property (3.9) is a consequence of the following identity, which is a 
special case of "Cartan's magic formula" for Lie derivatives and for a flow field 
generated by the dilations with center a. 

-^(tV(a + t(x - a)) = 

di-i(t^~^{x - a) J u(a + t{x - a)) j + - a) j u(a + t{x - a)) (4.6) 

Here u is an £-form. The result is 

d£_i R^u + R^+i d£U = u il<i<d-l); 

Ridou = u- (^,u)p^^ (^ = 0); (4.7) 

dd-i RdU = u {£ = d) . 

These relations are valid for all u G C^(M'^, A^) and by extension for all u E 

H'{D,A^), s e M. 

The perfect match of (4.7) with (3.9) from Assumption 4 suggests that the 
regularized Poincare lifting R^ provides suitable local liftings as stipulated in As- 
sumption 4. To this end, we can choose as local spaces of "more regular forms" 

X{K, A') := H{de, K) n H^{K, A') , 
S{K,A^-^):=H^{dt_i,K) and S{K,A^):=W\di,K), 

for some < r < 1, where we denote by H^{dk, K) the space 

W\dk.K) := {v e H'^{K,A'') : dfcv e H"^ {K , A'^+')} . 

All these spaces are equipped with the natural Hilbert space norms. Also keep in 
mind that the global spaces X{m, A^), S{Tl, A^^^) and S{Tl, A^) are determined 
by their local definition on the mesh cells K, cf. (3.3) and (3.6). For the particular 
choice (4.8) an assumption of Corollary 3.3 can be verified. 

Lemma 4.2 For X{m, A^) arising from (4.8) the space A^) n H^d^ 0, VL) 

o 

is dense in Hide, 0, u). 
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Proof. By [26, Thm. 4.9(c)] we have a direct decomposition 

if(d,0,fi) = d,„i H\n,K'-^) © , C C^(M", AO , (4.9) 

where C^(M'^, A^) is the space of compactly supported, smooth forms on 
with support contained in VL or, equivalently, the space of all smooth £-forms on Vt 
that vanish on dVt together with all their derivatives. Since C^(]R'^, A^~^) is dense 

in H^{VL, A^~^), we deduce: 

C^im.\K^)ndi^iH\VL,K^-^) is dense in d^^i ^^(fi, A^^^ 

As every u G C^(M'^, A^ belongs to X{M, A^, the assertion follows. □ 
We point out that the choice of r in (4.8) is determined by Assumption 1. 

Also note that whenever we opt for (4.8), Rellich's theorem ensures Assumption 

2, because the mesh is kept fixed. 

The construction of entails a constraint on the cell shapes. This is satisfied 

for standard finite element meshes, where the cells usually are convex polyhedra. 

Assumption 9 Every cell G SOT is a star-shaped polyhedron. 



Lemma 4.3 Assumption 9, the choice (4.8)/or spaces X{K, A^) and S{K, A^ ^ 
imply Assumptions 2, 3 and 4. 

Proof. The only fact which remains to be proved is the maximality relation in 
Assumption 3 

d^_i S{m, A^-^ = d^_i H(de_i, Q) n x{m, A^. 

The inclusion C holds by definition. Let us prove the converse inclusion. 
Let u G d^_i H{de^i, fi) n X(9Jt, A^. Thus u = d^„i </> with G H{de^i, fi). 
Since u G A^, using [26, Cor. 4.7] we obtain that there exists ■j/' e H^{Q, A^^^ 

such that u = d£„i -j/j. In particular, ?/? belongs to H'^{K, A^~^) for all K and, 
since u |^^, belongs to H''{K, A^, we finally find that i^l^^e H''{de^i, K). □ 

4.4 Preservation of polynomial forms 

Fundamental in finite element methods is the notion of polynomial differential 
forms. For an ordered £-tuple / = (zi, . . . , ig), ii < 12 < ■ ■ ■ < in, ■ ■ ■ C 
{l,...,4,let 

dxi := dxi^ A ■ ■ ■ A dxi^. , 
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where dxj, j = 1, . . . , d, are the co-ordinate 1 -forms in Euclidean space M''. The 
space 'Pp(A^) of polynomial ^-forms on R'^ is defined as 

Vp{A') := {u = ■ e Pp(M'')} , 

where indicates summation over all ordered ^-tuples, and Vp(W^') is the space 
of (i-variate polynomials of total degree < p. We remark that for d E {2, 3} 
polynomial forms possess polynomial vector proxies. 

From the definition (4.2) it is clear that the Poincare operator „ maps dif- 
ferential forms with polynomial coefficients to differential forms with polynomial 
coefficients. The same holds for the regularized Poincare operator by (4.3). 
If we want R^ to map a space P(A^) of differential forms of order £ (e.g., with 
polynomial coefficients) into a space P(A^^^) of differential forms of order £ — 1, 
it is sufficient to require the following two properties, see [26, Proposition 4.2]. 

Proposition 4.4 Assume that P(A^) and P(A^~^) are finite- dimensional spaces 
of differential forms satisfying 

(i) The space P(A^) is invariant with respect to dilations and translations, that is 

For any teR,aeR''' : ifu G P(A^), then {x ^ u{tx + a)) G P(A^) . 

(ii) The interior product xj : u x ju maps P(A^) to P(A^~^). 
Then Re maps P(A^) into P(A^~i). 

For the compatibility Assumption 5 to hold, it is therefore sufficient to make 
the following assumption about the local polynomial space Vp{K). 

Assumption 10 

(i) The space Vp{K) is invariant with respect to dilations and translations. 

(ii) The differential operator a;jd^:ui— T-xjd^u maps Vp{K) into Vp{K). 

To summarize: 

Assumptions 9, 10, and (4.8) =^ Assumptions 2, 3, 4, and 5 

5 Discrete differential forms 

Now we introduce concrete spaces of discrete differential forms. We merely sum- 
marize the constructions that have emerged from research in differential geometry 
(the "Whitney-forms" introduced in [48]) and finite element theory ("Raviart- 
Thomas elements" of [47] and "Nedelec finite elements" of [44, 45]). These 
schemes were later combined into the concept of discrete differential forms [14, 
36]. Surveys and many more details can be found in [37, 4, 5, 15]. 
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5.1 Simplicial meshes 



Let SOT be a conforming simplicial finite element mesh covering the Lipschitz 
polyhedron c W^. As elaborated in [4, Sect. 3 & 4] for p G N the following 
choices 

V'^{K) ■.= V,^i{A')\^ + x^Vp^i{A'+')\^ (5.1) 

and 

V^iK) := P,(AO 1^ (5.2) 
of local spaces, through (3.1), gives rise to meaningful global finite elment spaces 

O 

V^IJXH) of discrete differential forms. 

By construction both Assumption 9 and Assumption 10 are satisfied for these 
spaces. The asymptotic density property also holds. 

Lemma 5.1 The spaces V^{Vyt) of discrete differential forms built from (5.1) or 
(5.2) meet the requirement (2.12). 

Proof. It is a classical result of finite element theory that the spaces of degree 
p Lagrangian finite element functions Vp(97t) are asymptotically dense in H^iVt). 
Thus the space of polynomials £-forms with coefficients in Vp(OJt), which is a 
subspace of Vp(5[Jt), is asymptotically dense in i/^(f2,A^). The latter space is 

O . 

obviously dense in //(d^, Q), since this is already true for C^(f2, F'). □ 



5.2 Tensor product meshes 

Let SOT be a conforming finite element mesh of the Lipschitz polyhedron ^2 whose 
cells are affine images of the unit hypercube K in W^: for ii' G SOT the we write 
#x : K ^ K for the associated unique affine mapping. We generalize the 
construction of [44]: on the cube we define (with notations introduced in Section 
4.4) 

Vp(^) := |9 = ^uidxi, ui{x) = Yluijixj), uij G j^''"^ ^ / } ' 

The local spaces are obtained by affine puUback 

V^W := {^K-)X(^) ■ (5.3) 

This affine tensor product construction also complies with Assumption 9 and As- 
sumption 10. Completely parallel to Lemma 5.1, one proves the following result. 
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Lemma 5.2 The requirement (2.12) is satisfied for the spaces V^(OJt) spawned 
by (5.3). 

Remark 5.3 For all the above meshes the cells are affine images of a single ref- 
erence cell, the "unit simplex" or "unit hypercube". We could allow non-affine 
cells, because, under the assumption that the transformations are "nearly affine", 
see [22, §4.3], all crucial estimates like Lemma 3.1 can be transferred to the ref- 
erence cell using the puUback of differential forms. Of course, also the projection 
operators vTp would have to be defined by transformation. □ 

6 Application in dimensions two and three 

We adopt the discrete spaces from Sect. 5 along with the regularized Poincare 
lifting from Sect. 4. We rely on the choice (4.8) for spaces X and S, with a 
regularity exponent re (0,1] which has to be chosen suitably. 

In order to establish the discrete compactness property from Definition 2.3, it 
remains to verify the regularity Assumption 1 and Assumptions 6, 7, and 8 for 
convenient local projectors vr^ ^. 

Local projectors which make the discrete diagram of Assumption 7 commute 
do exist in the general framework of differential forms of any degree. They gen- 
eralize Nedelec edge element projections and can be referred to as moment based 
projection operators. They are suitable for the /i-version of finite elements in di- 
mensions 2 and 3. In higher dimensions some of them (for low degree forms) 
require a higher regularity than H'^ to be defined. In [4, 5], they are modified 
by an extension-regularization procedure in order to be defined and bounded on 
L^. However, such operators cannot be used for the p version of finite elements, 
because no estimates (stability or error bounds) are known with respect to the 
polynomial degree p. 

The proper projection operators for p- version approximation are so-called pro- 
jection based interpolation operators, see [28, 29, 18, 27, 30]. Variants for any i 
and d are available and they are designed to commute in the sense of Assumption 

7 [37, Sect. 3.5]. 

At this point we have to abandon the framework of general i and d, because 
both regularity results (Assumption 1) and the analysis of projection operators 
(Assumption 8) are not presently available for general i and d. We have to dis- 
cuss them for special choices of i and d separately, relying on a wide array of 
sophisticated results from the literature. 

Theorem 6.1 (Convergence of Galerkin approximations) For d = 2, 3, and < 
i < d, the Galerkin discretization of (2.5) on a Lipschitz polyhedron based on any 
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of the families of discrete differential forms introduced in Sect. 5 ojfers a spectrally 
correct, spurious-free approximation. 

Proof. We skip the case £ = 0, for which the standard Galerkin approximation 
theory for operators with compact resolvent can be applied, see [41]. 

To begin with, we focus on the discrete compactness property and verify the 
assumptions 1, 6, 7, and 8 for d = 2 and = 3 separately. 

• d = 2, (. = 1: m terms of vector proxies we find the correspondence 

F(di,fi) ~ i?(curl,l]) nif(divO,fi) . (6.1) 

Regularity theorems for boundary value problems for — A on the polygon confirm 
the existence of 5 = 5{Vt) > such that 

H{cur\, n) n H{d\y 0, fi) C H^+^/^{n) , (6.2) 

in the sense of continuous embedding, see [34, Sect. 3.2]. This suggests to choose 
r = (5+1/2 in (4.8) and Assumption 1 will hold true. Hence, we deal with the 
concrete spaces 



X{m, A^) = H{cur\, n)nYl {H^^^'^{K)f , (6.3) 

Kern 

5(9Jl,A°) = H\n) n Y[ H^+^/\K) . (6.4) 

Commuting local projection based interpolation operators tt^ and vr^ have 
been proposed for triangles and for quadrilaterals in [28]. With the choice (6.3) 
and (6.4) they live up to Assumptions 6 and 7. Assumption 8 holds with Eo{p) = 
Cp~^/'^ and C > depending only on the shape-regularity of the cells, cf. [27, 
Thm. 4.3] and [6, Thm. 4.1]. Finally, these interpolation operators satisfy the nat- 
ural condition of conformity (3.14) by construction, which makes they meet all 
our requirements, c/ Sect. 3.4. 

• (i = 3,£ = l,2:we have the vector proxy incarnation 

yid.n)-/^''""-"'^^*^'""-"' f"^ = ^- (6.5) 

^ ' |/f(div,fi) nif(curlO,fi) for£ = 2. 

Citing results from [1] and [37, Sect. 4.1], we find 5 = 5{^l) g]0, |] and continu- 
ous embeddings 

^(curl, n) n if(div, fi), ^(div, fi) n if(curl, il) c H^^^^'^{n) . (6.6) 



27 



Therefore, using the construction (4.8) with r = 5+1/2, Assumption 1 is satisfied 
for i E {1,2}. The relevant spaces of more regular forms now read 





= H{cm\,Q) n Yl {H^+^'^{K)f , 
Kern 


(6.7) 


X(9Jt,A') 


= i?(div,fi) n JJ {H^+^'^{K)f , 


(6.8) 


s{m, A°) 


Kern 


(6.9) 




= H{cur\,Q) n Yl H^+^/^{cur\,K) . 


(6.10) 



The essential commuting local projection based interpolation operators tt^^, 
m = 0, 1, 2, have been introduced in [29] for tetrahedral meshes and in [28] for 
meshes comprising parallelepipeds. By construction they comply with Assump- 
tion 7. Assumption 6 for the spaces S{M, A°) and S{Tl, A^) from (6.9) and 
(6.10), respectively, and r = 5+ l/2isa consequence of Sobolev embedding the- 
orems. Relying on [27, Th.5.3] we obtain like in the 2D case that in Lemma 3.1 
we can take Emip) = Cp~^^'^ for m = and m = 1. 

• Finally, we appeal to Lemmas 5.1, 5.2 together with Lemma 4.2 and apply the 
abstract theory of Sect. 3 in the form of Corollary 3.3 to conclude the proof of the 
theorem. □ 

Corollary 6.2 (Approximation of the Maxwell eigenvalue problem) The p ver- 
sion finite element discretization of the Maxwell eigenvalue problem (1.2) based 
on edge elements from the first or second Nedelec family on triangles or on tetra- 
hedra, or from the first Nedelec family on parallelograms or on parallelepipeds 
offers a spectrally correct, spurious-free approximation. 

Remark 6.3 Instead of (1.2) we may consider the variational formulation of the 
more general Maxwell eigenvalue problem (1.1), corresponding to the case of 
anisotropic inhomogeneous material: 

Seek u G ^(curl, Vl) \ |0|, G Mn such that 

(6 11) 

[jj,'^ curl u, curl v)^2(j^) = ^'^ (eu, v)^2(f^) Vv G i3"(curl, i7) , 

with uniformly positive material tensors fi = /i(x), e = e(x). The same edge 
element discretizations listed in Corollary 6.2 provide spectrally correct, spurious- 
free approximations of this problem. This generalization of Corollary 6.2 can be 
achieved with standard tools (see, in particular. Propositions 2.25, 2.26, and 2.27 
of [19], and [38, Sect. 6], [37, Thm. 4.9]). □ 
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Remark 6.4 The restriction on the families of elements mentioned in the Corol- 
lary is essentially due to the availability of published results about suitable inter- 
polation operators. Thus, for example, as soon as a generalization of the p version 
error estimates of [29, 27] for projection-based interpolants to meshes containing 
prismatic or more general polyhedral elements becomes available, our result about 
the approximation of the Maxwell eigenvalue problem will apply to such meshes, 
too. □ 

Remark 6.5 Several obstacles prevent us from establishing the assumptions of 
the abstract theory for > 3. On the one hand, continuity properties of projection 
based interpolation operators have not been investigated for d > 3. Also, regu- 
larity results along the lines of (6.6) are have not been published for polyhedra in 
higher dimensions. 

On the other hand, the innocuously looking requirement (3.14) for the projection 
operators — corresponding to the requirement that the global projection operators 
are constructed elementwise from local degrees of freedom — entails that the trace 
of forms in S{K, A^~^) onto ^— 1-dimensional facets in 5^£„i(A') must make sense. 
However, we cannot expect more than regularity for the space S'(A', A^~^). 
Hence, by trace theorems for Sobolev spaces, the spaces ^(A', A^~^) allow for 
traces on m-facets for m > | — 2 at best, which means that £ > | — 1 is required 
to allow for the construction of a local projection based interpolation complying 
with Assumptions 6 and 8. 

Perhaps, an analysis in L^-spaces as in [1, Lemma 4.7] can make possible an 
extension of the theory to higher dimensions, but this is beyond the scope of the 
present article. □ 

Remark 6.6 Our approach does not cover /ip-refinement, for various reasons. 
One reason is that there exist many variants of hp refinements in 3D, and covering 
them would in any case require a much longer paper than the present one. 
Another reason is technical: The existing convergence proof of the hp approxima- 
tion of the Maxwell eigenvalue problem in [11] — while also based on the proof 
of the discrete compactness property — uses a different technical tool, namely an 
estimate of the stability of a certain projection operator. This kind of estimate 
is currently only available for intervals in ID and for rectangles in 2D. 
The technique used in the present paper is based on the regularized Poincare lift- 
ing, and adjusting this to variable polynomial degree poses formidable technical 
challenges. Only in 2D these could be mastered so far, as was demonstrated in [7] 
in the context of boundary element analysis. □ 
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7 Conclusion 



In this paper we have proved that the p-version of finite elements based on gener- 
alized Nedelec edge elements provides a spurious-free spectrally correct approx- 
imation of the Maxwell eigenvalue problem. The essential point was the proof 
of the discrete compactness property. We showed that this property follows from 
a set of rather natural assumptions about the family of finite element spaces and 
interpolation operators, and in addition we showed that these assumptions are 
implied by recently found results on lifting operators and on projection-based in- 
terpolants. 

In the approach pursued in [4, 5] the discrete compactness property is not 
addressed directly: in the framework of the /i-version for differential forms, mod- 
ified moment-based projection operators are used. These new operators satisfy the 
strong property of being uniformly bounded in and are constructed by means 
of a delicate extension-regularization procedure, see also [21, 20]. 

On the one hand this uniform boundedness property is stronger than our as- 
sumption 6 and replaces in a certain way the discrete compactness property. But 
on the other hand, it is currently not known whether a construction of projec- 
tion operators by extension-regularization could also be employed in the case of 
the p-version of finite elements, or whether the construction of a p-uniformly L^- 
bounded family of cochain projections is even possible. 
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